Mechanism for spectral break in cosmic ray proton spectrum from Supernova remnant W44 
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Recent observations of the supernova remnant W44 by the Fermi spacecraft observatory strongly support the 
idea that the bulk of galactic cosmic rays is accelerated in such remnants by a Fermi mechanism, also known 
as diffusive shock acceleration. However, the W44 expands into weakly ionized dense gas, and so a significant 
revision of the mechanism is required. In this paper we provide the necessary modifications and demonstrate 
that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of 
accelerated particles by exactly one power. The spectral break is caused by Alfven wave evanescence leading 
to the fractional particle losses. The gamma-ray spectrum generated in collisions of the accelerated protons 
with the ambient gas is also calculated and successfully fitted to the Fermi Observatory data. The parent proton 
spectrum is best represented by a classical test particle power law oc E^^, steepening to E^^ at E/j^ ~ IGeV due 
to deteriorated particle confinement. 



Ongoing direct observations of a number of galactic super- 
nova remnants (SNRs) in the X- and gamma-ray bands sup- 
ported by the radio, optical and UV data fT,*?] seem to rapidly 
close in on the century long problem of the cosmic ray (CR) 
origin. Overall, the observations favor the diffusive shock ac- 
celeration (DSA, operating in SNR shocks) as a mechanism 
for the production of galactic CRs. However, there are at 
least two critical questions that observations pose to the the- 
ory. First, what exactly is the form of the spectrum that the 
theory predicts? Second, when do we really see the proton 
(i.e., the primary CR component) emission and when do we 
see just a contaminating but radiatively more efficient elec- 
tron component? 

Several decades of studies of the background galactic CRs, 
preceding direct observations (that became available only re- 
cently) demonstrated that already a simple test particle ver- 
sion of the DSA (leading to a power-law ^ E^^ particle en- 
ergy distribution) reproduces the form of the energy spectrum 
reasonably well. It is possible, however, that a somewhat 
steeper power-law would better accommodate the E^^-'' mea- 
sured spectrum [HI 01 (see also the Discussion section below), 
even though this difference can be attributed to CR propaga- 
tion losses (also not completely certain). Nevertheless, the 
direct observations of emission coming from particles accel- 
erated in SNR shocks often indicate significantly steeper than 
E^^ spectra [SJ. Also the recent Fermi Large Area Telescope 
(LAT) observations of the SNR W44 suggest that the spectrum 
of the gamma ray producing protons is substantially steeper at 
high energies than the DSA predicts |6]. It should be noted 
that a similar discrepancy has already been found in the high 
energy gamma ray spectrum of the SNR RXJ 1713 measured 
by the CANGAROO |7] and HESS |5] atmospheric Cerenkov 
telescopes. The simple DSA prediction was even used as an 
argument against the hadronic origin of the gamma emission 
on the ground that if the high energy part of the spectrum was 
extrapolated (with the DSA index 2) to lower energies, the 
emission there would be unacceptably high, e.g. [8J. Moti- 
vated by this inconsistency with the DSA theory in particular, 
we have previously published a suggestion of the mechanism 
for the spectral break |9], which can resolve the discrepancies. 
We believe that a similar phenomenon is also applicable to the 
recent observations of the SNR W44. 



In particular, we show below that ion-neutral collisions in 
a molecular cloud adjacent to the remnant steepen the energy 
spectrum of accelerated particles. We calculate the gamma- 
ray spectrum resulting from the collisions of accelerated pro- 
tons with the molecular cloud material and fit the result to the 
Fermi Observatory data. 



Results 

Mechanism for the spectral break The physics of the 
spectral break is very simple. When a SNR shock approaches 
a molecular cloud (MC) or a pre-supernova swept-up shell, 
confinement of accelerated particles generally deteriorates. 
Due to the particle interaction with magnetic turbulence, con- 
finement requires scales similar to the particle gyroradius 
ifioi [Till . While the waves are in a strongly ionized (closer 
to the shock) medium they propagate freely in a broad fre- 
quency range at the Alfven speed Va = B/ \/4npj with the 
frequencies co — Wa- Here k is the wave number (assumed 
parallel to the local field B) and p, is the ion mass density. 
As long as the Alfven wave frequency is higher than the ion- 
neutral collision frequency v,„, the waves are weakly damped. 
When, on the other hand, the ion-neutral collision frequency 
is higher (deeper into the cloud), neutrals are entrained by the 
oscillating plasma and the Alfven waves are also able to prop- 
agate, albeit with a factor y/pi/po < 1 lower speed, where 
Po is the neutral density. The propagation speed reduction 
occurs because every ion is now "loaded" with po/p, neu- 
trals. In between these two regimes Alfven waves are heavily 
damped and even disappear altogether for sufficiently small 
Pi/Po ^0.1. The evanescence wave number range is then 
bounded by ki — v,„/2Va 



and k2 = l^/pJpnVm/VA- These 
phenomena have been studied in detail in |12L[l3n . and specif- 
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ically in the context of the DSA in [14-16]. Now we turn to 
their impact on the particle confinement and emissivity. 

In the frame work of a quasilinear wave-particle interaction 
the wave number k is approximately related to the parallel (to 
the magnetic field) component of the particle momentum 
by the cyclotron resonance condition kp\\ jm = ±(0c where the 
(non-relativistic) gyro-frequency (Oc = eB/mc. Therefore, the 
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frequency range where the waves cannot propagate may be 
conveniently translated into the parallel momentum range 

Pi < \P\\ \ < P2 (1) 

with 

pi = IVAinCOc/Vin , f 2 = ^ VPo/Pi > PU (2) 

That a spectral break must form at the photon energy corre- 
sponding to the particle momentum p~p\— Pbn can be read- 
ily understood from FiglT] The 'dead zones' pi < |/?|| | < p2 
imply that particles with | > pi do not turn around (while 
moving along the magnetic field) and escape from the region 
of CR-dense gas collisions at a p\\/p fraction of the speed of 
light. This happens in the region of enhanced gamma radia- 
tion. Therefore, an initially isotropic distribution of acceler- 
ated particles fills only a slab in momentum space |/?|| | < pi 
and becomes highly anisotropic (a 'pancake' distribution). 
What matters for the integral emission, however, is a locally 
isotropic component / of this new proton distribution. It can 
be introduced by re-averaging the 'pancake' | < pi) dis- 
tribution in pitch angle, f (p) = Jq f{p,^)d^, and is read- 
ily obtained assuming that particles remaining in the dense 
gas (those with | < pi) maintain their isotropic pitch-angle 
distribution, i.e. 

Ml 



where /o (p) is the initial (isotropic) distribution function 
and ji is the cosine of the pitch angle (see FiglU, jUi — 
min{;?i /p, 1}. Thus, the slope of the particle momentum dis- 
tribution becomes steeper by exactly one power above p = 
pi = phr- In particular, any power-law distribution p^'', 
upon entering an MC, turns into p^''^' 3t p > Pbr^ and pre- 
serves its form at p < phr- 

Now the question is whether particle escape from the MC 
can quench the acceleration process itself. In principle it can, 
since also the particles from the interval | < pi cannot re- 
turn to the shock and continue acceleration. Instead, by slowly 
leaking to the part of the momentum space with | > pi, 
they should ultimately escape as well. This would certainly 
be the case if the MC were filling the entire shock precursor 
However, MCs are known to be clumpy iflTl - fTgll . and may fill 
only a small fraction of the precursor In this case the accel- 
eration process continues almost unimpeded, the accelerated 
protons illuminate the 'cloudlets' and make them visible in 
7-rays due to the high density target material. 

Break momentum While the one power spectral break in 
the pitch-angle averaged particle distribution seems to be a 
robust environmental signature of a weakly ionized medium 
into which the accelerated particles propagate, the break mo- 
mentum remains uncertain. According to eq.dU, pi,r (= p\) 
depends on the magnetic field strength and ion density as well 



as on the frequency of ion-neutral collisions, v,„ = no (o'V). 
Here {oV) is the product of the collision cross-section and 
collision velocity averaged over the thermal distribution. Us- 
ing an approximation of 1151 l20ll for (crV), pf,, can be esti- 
mated as 

Phr/mc ~ 10b2 r^-O-^HQ i„-^/2^ (4) 

Here the gas temperature is measured in the units of IQ'^K, 
magnetic field B^x -in microgauss, uq and n,- (number densi- 
ties corresponding to the neutral/ion mass densities po and 
Pi) -in cm^^. Note that the numerical coefficient in the last 
expression may vary depending on the average ion and neu- 
tral masses and can be higher by a factor of a few for typical 
molecular cloud conditions lfT2il2Tll than the estimate in eq.(|4]l 
suggests. The remaining quantities in the last formula are also 
known too poorly to make an accurate independent prediction 
of the position of the break in the gamma ray emission region. 
Those are the regions near the blast wave where complicated 
physical processes unfold. They include particle acceleration, 
strong MHD turbulence (driven by particles and their inter- 
action with ambient gas inhomogeneities), gas ionization by 
shock generated UV photons, turbulent plasma heating and 
even evaporation of magnetic cloudlets |20, 22l.|23[]. Also im- 
portant may be the ionization by the low energy CRs acceler- 
ated at the blast wave. However, as their diffusion length is 
shorter than that of the particles with p > pt,r, we may assume 
that they do not reach the MC. Pre-ionization by the UV pho- 
tons can also be ignored for the column density > 10 ^^cm^^ 
ahead of the shock beyond which they are absorbed 02411 . The 
authors or the Ref. f24], using the earlier data from f25'l have 
also analyzed the parameters involved in eq.© and found the 
above estimate of pi„ to be in a good agreement with the spec- 
tral break position measured by the Fermi LAT. Nevertheless, 
we may run the argument in reverse and use the Fermi obser- 
vations 1 6] of the gamma-ray spectrum of SNR W44 to de- 
termine the break momentum in the parent particle spectrum 
and constrain the parameters in eq.®. Once we also know 
the amount of the slope variation Aq, we can calculate the full 
spectrum up to the cut-off energy. 

Particle spectra To calculate the particle spectra, we need 
to determine the degree of nonlinear modification of the shock 
structure. In principle, it can be calculated consistently, given 
the shock parameters and the particle maximum momentum, 
Pmax- In the case of a broken spectrum, pi,r plays the role of 
Pmax, since this is the momentum where the dominant con- 
tribution to the pressure of accelerated particles comes from, 
thus setting the scale of the modified shock precursor. Note 
that in the conventional nonlinear (NL) acceleration theory, 
the cut-off momentum p,nax plays this role, because the non- 
linear spectra are sufficiently flat so as to make the pressure 
diverge with momentum. 

The break in the photon spectrum is observed at about 2 
GeV, which places the break on the proton distribution at 
about Phr — IGeV/c Ig]. For the strength of the break A<7 = 1, 
the spectrum above it is clearly pressure converging, and we 
perform the calculation of the shock structure and the spec- 
trum using this break momentum as the point of the maximum 
in the CR partial pressure. Once this point is set, we can use 



3 



an analytic approach for a stationary nonlinear acceleration 
problem using as an input parameter 

Apart from pj,^, the nonlinear solution depends on a num- 
ber of other parameters, such as the injection rate of thermal 
particles into acceleration, Mach number, the precursor heat- 
ing rate and the shock velocity V,. Of these parameters only 
the latter is known accurately, Vs ~ 3QQkm/s, the other param- 
eters are still difficult to ascertain. Fortunately, in sufficiently 
strong shocks the solution generally tends to either stay close 
to the test particle (TP) solution (leaving the shock structure 
only weakly modified) or else it transitions to a strongly mod- 
ified NL-solution regime. The TP regime typically establishes 
in cases of moderate Mach numbers, low injection rates and 
low p„,ax (now pi,r), while the NL regime is unavoidable in 
the opposite part of the parameter space. 

In the TP regime the spectrum is close to a power-law 
with the spectral index 2 throughout the supra-thermal energy 
range. In the NL regime, however, the spectrum develops a 
concave form, starting from a softer spectrum at the injection 
energy, with the index ^ ~ {rs+2)/{rs— 1 ) > 2, where < 4 
is the sub-shock compression ratio. Then it hardens, primarily 
in the region p ^ mc, where both the partial pressure and dif- 
fusivity of protons change their momentum scaling. The slope 
reaches its minimum at the cut-off (break) energy, which, de- 
pending on the degree of nonlinearity, can be as low as 1.5 
or even somewhat lower if the cut-off is abrupt. The question 
now is into which of these two categories the W44 spectrum 
falls? Generally, in cases of low maximum (or, equivalently, 
low spectral break /^^^ ^5 10) momentum, the shock modifi- 
cation is weak, so the spectrum is more likely to be in the 
only slightly nonlinear, almost TP regime. On the other hand, 
there is a putative indication from the electron radio emission 
that their spectrum may be close to qe ~ 1.75, which could 
be the signature of a moderately nonlinear acceleration pro- 
cess. It should be remembered, however, that this is a global 
index across the W44 remnant. There are resolved bright fil- 
aments where a canonical a = —0.5 spectrum, corresponding 
precisely to the TP parent electron spectrum with = 2 is 
observed i26l l. Moreover, there are regions with the positive 
indices a < 0.4 which cannot be indicative of a DSA pro- 
cess without corrections for subsequent spectral transforma- 
tions such as an absorption by thermal electrons |26]. These 
regions can very well contribute to the overall spectral harden- 
ing discussed above, mimicking the acceleration nonlinearity. 
Finally, secondary electrons give rise to the flattening of the 
radio spectrum as well Il24ll . 

If the accelerated protons and electrons respond to the 
turbulence similarly, which is almost certainly the case in 
the ultra-relativistic regime, their spectra should have similar 
slopes there (as long as the synchrotron losses are ignorable). 
In using the electron radio spectrum as a probe for the level of 
acceleration nonlinearity, the following two relations are use- 
ful. First, there is a relation between the electron energy and 
the radio frequency Vmhz = 4-6 • B^E^^y. The second rela- 
tion, ^f, = 1 — 2a, links the spectral index of radio emission a 
(assuming the radio flux v") and the spectral index of the 
parent electrons qe (assuming their energy spectrum oc E^^ll- 
Once the global radio spectral index of W44, a ~ —0.37 12611 



is generated by freshly accelerated electrons in the frequency 
range 74 < v < 10700 MHz, the electrons should maintain 
their modified spectrum over the energy range spanning more 
than one order of magnitude. For example, assuming ~ 70 
[6], one sees that electrons must maintain an index q^ « 1.75 
between 0.46 < E < 5.8 GeV. While the upper bound is ac- 
ceptable given the spectral break proton energy inferred from 
the super GeV emission measured by the Fermi LAT, the 
lower end is rather uncomfortable, since the nonlinear modifi- 
cation of both protons and electrons with the Bohm (or other 
similar for protons and electrons turbulent diffusivities) starts 
(slowly) only at the proton rest energy. The calculated non- 
linear spectra are shown in Fig |2] for the both species. At and 
below IGeV, the electron spectrum is very close to the test par- 
ticle solution, qe ^ 2, even though the proton spectrum may 
be somewhat steeper there, as we mentioned above. 

Photon spectra The above considerations somewhat 
weaken the radio data as a probe for the slope of the elec- 
tron and (more importantly) for the proton spectrum. There- 
fore, the exact degree of nonlinearity of the acceleration re- 
mains unknown and we consider both the TP and weakly NL 
regimes in our calculations of the photon spectra, generated 
in p — p collisions. The results are shown in Fig|3] The best 
fit to the Fermi data is provided by a TP energy distribution 
(oc E^^) below phr — 7GeV/c with the spectrum steepening 
by exactly one power above it. Note that the small devia- 
tion of the computed spectrum from the lowest energy Fermi- 
point may also be corrected by including the Bremsstrahlung 
of the secondary electrons [24]. The spectrum steepening is 
perfectly consistent with the proton partial escape described 
above (with no parameters involved) and shown in FiglT] For 
comparison, a weakly NL spectrum, shown in Fig|2l is also 
used for these calculations (dashed line in FigO, but its fit 
would require a stronger break (A^ > 1) or a low momentum 
cut-off, Fig[3] i.e. at least one additional free parameter 



Discussion 

To summarize the results, the previously suggested mecha- 
nism for a break on the spectrum of shock accelerated protons 
|9] is found consistent with the recent |6] Fermi LAT obser- 
vations of the SNR W44. The observed gamma ray spectrum 
most likely results from the decay of ;r"-mesons which are 
born in p — p collisions of shock accelerated protons with an 
ambient dense gas. The parent proton spectrum is best rep- 
resented by a classical test particle power law E^^, steep- 
ening to E^^ at Ehr ~ IGeV due to deteriorated particle con- 
finement caused by the ion-neutral collisions and the resultant 
Alfven wave evanescence. The position of the break momen- 
tum in the particle spectrum may be estimated using eq.©, 
or conversely, the combination of parameters involved in this 
estimate can be inferred from the measured break momentum. 
The cut-off momentum is not constrained in this scenario. 

An alternative explanation, based on a different mechanism 
of the break, associated with the change of the particle trans- 
port in the CR shock precursor |27] is also possible but is less 
definitive in the spectrum slope variation ls.q across the break 
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(see also |* 24|] for the most recent alternative suggestions). In 
addition, the ref fz?"] mechanism would imply a considerable 
nonlinearity, i.e. a stronger CR shock precursor than it seems 
appropriate for the observed low break momentum. Still alter- 
natively, assuming our "environmental" break mechanism is 
at work, i.e. Aq = 1, but the shock structure appreciably mod- 
ified, we arrive at the Zi^'-^^ spectrum below the break (as 
the radio observations suggest for the electrons), and Zi^^^^ 
above the break. A fit to the data is marginally possible, but it 
would require a relatively low cut-off momentum at about 100 
GeV/c. This possibility may be supported or ruled out once 
the data (upper limit) around this energy become available. 

The most robust and attractive aspect of the suggested 
mechanism for the spectral break is the exact Aq = 1 varia- 
tion of the spectral index. Indeed, this change in the spec- 
tral slope is due to truncation of the particle momentum phase 
space and does not depend on any parameters. In a combina- 
tion with the test particle regime operating below the break, 
which is consistent with the low value of the break momen- 
tum, it provides a very good fit to the Fermi LAT data. Of three 
physically different types of spectral breaks suggested earlier 
(9, 27, 28], namely the current, "environmental" mechanism 
appears to be plausible in situations where dense target gas 
is present, as required for efficient ;r*' production. In a more 
general context, a spectral break is a natural resolution of the 
well known but puzzling trend of the nonlinear (i.e. suppos- 
edly improved) DS A theory to develop spectra which are even 
harder than the simple test particle spectra, thus encountering 
more difficulties while accommodating the bulk of observa- 
tions |3, 4]. Such a spectrum - i.e., diverging in energy- ex- 
hausts the shock energy available for the acceleration as the 
cut-off momentum grows. Note, that the spectrum of the RX 
J 17 13.7-3946 [5] is also consistent with the same break mech- 
anism but with phr ^ IQ^GeV /c and with naturally stronger 
acceleration nonlinearity However, it is difficult to make 
the case for hadronic origin of the gamma-ray emission of the 
RX J1713.7-3946 S H S. The fundamental role of the 
W44 for the problem of CR origin is that this particular rem- 
nant seems to rule out contaminating electron emission due 
to Bremsstrahlung and the inverse Compton scattering ll6l l24ll 
thus favoring the hadronic origin of the gamma emission. 

Methods 

Acceleration model and proton spectrum Methods of cal- 
culation of the accelerated particle spectra, including the par- 



ticle back-reaction on the shock structure, are now well devel- 
oped. We use the diffusion-convection equation for describ- 
ing the distribution of high energy particles (CRs). To include 
the back-reaction, three further relations are used in a quasi- 
stationary acceleration regime. These are the conservation of 
mass and momentum fluxes in the smooth part of the shock 
transition (so called CR-precursor) and the Rankine-Hugoniot 
relation for the shock compression at the gaseous disconti- 
nuity (sub-shock). The complete system of equations is then 
reduced to one nonlinear integral equation 1 3 1 ] . This equation 
is solved numerically and the results are shown in Fig|2] The 
above method of calculation of the parent particle spectra is 
shown to be in excellent agreement |32] with numerical simu- 
lations ll33ll34ll as well as with other semi-analytic approaches 

Gamma-ray emission spectrum Once the parent proton 
spectrum is obtained, we calculate the ;r*' production rate and 
the gamma-ray emissivity. In so doing, we adopt numerical 
recipe described in detail in f3^,'37]. The physical processes 
behind these calculations are [38] (i) collisions of accelerated 
protons with the protons of the ambient gas resulting in the 
following spectrum of tt" -mesons: 



F,, = AKNp, I ^^iEllEAj^^ (Ep) dEp (5) 



where, Npf, is the number density of protons in the gas, 
da /dEji is the differential cross section for the production 
in collisions between accelerated protons of energy Ep and 
gas protons, Jp is the flux of accelerated protons, is the 
energy of Tt" mesons; (ii) decay of Ti^ resulting in the gamma 
emission spectrum 



Ey+mlc'^/4Ey 



where m;^ is the pion rest mass. 
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Figure 1: Momentum space of accelerated protons. Particle scat- 
tering zones on the (j>\\,p±^ - plane of momentum space. Pr in the 

stripes p\ < I < p2 are not scattered by waves (see text). There- 
fore, particles from the domains | | > P2 maintain their propagation 
direction and promptly escape from the dense gas region. 
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Figure 2: Spectra of accelerated protons and electrons. The both 
particle distributions are calculated for a weakly modified shock 
and are shown in momentum normalization (f (p) is steeper by 
two powers than the spectra in energy normalization, used in the 
text). Both spectra are multiplied by p^, so that the test parti- 
cle distribution is flat. Shock parameters: acoustic Mach num- 
ber M = 30, shock velocity Vs/c = 10^"', the break momentum 
Ph,. ~ line. Shock pre-compression (flow compression across the 
CR precursor) R=1.8, injection parameter v ~ 0.1 [defined as v = 
(4;r/3) [mc^/pV}) [pinj/mc]^ f {pi„j) , with p and Vj being the am- 
bient gas density and the shock speed, respectively]; injection mo- 
mentum pinj/mc ~ 1.4- 10^"'. 
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Figure 3: Gamma radiation spectra. Photon spectra resulting from 
decay and calculated for two different parent proton spectra com- 
pared against the Fermi data (open circles). Solid line: a test particle 
acceleration regime with the spectral index q = 2 below the break 
and (? = 3 above the break at pj,r = IGeVlc. Dashed line: a mod- 
erately nonlinear acceleration regime corresponding to the spectrum 
shown in Fig|2](i; ~ 1.75 and q ~ 2.75 below and above the break, 
respectively). Cut-offs are placed at 300 GeV for TP- and 100 GeV, 
for NL-spectrum. 



